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SECTION  I 


INTRODUCTION  AND  SUMMARY 


This  report  documents  the  results  of  investigation 
of  the  nuclear  weapons  effects  module  of  TACWAR,  a theater- 
level  computerized  warfare  simulation  code*,  with  the  pur- 
pose of  verifying  the  accuracy  of  its  data.  Periodic 
checks  of  nuclear  effects  data  in  wargaming  codes  are 
desirable  to  ensure  that  errors  do  not  cause  discrepancies 
in  the  intercomparison  of  various  models  and  to  preclude 
the  propagation  of  erroneous  results  or  the  influencing  of 
decisions  affecting  force  posture  or  hardware  procurement 
based  on  faulty  nuclear  effects  data. 

This  report  is  not  intended  to  stand  apart  from 
the  extensive  body  of  TACWAR  documentation.  It  assumes 
reader  familiarity  with  the  capabilities  and  content  of 
the  TACWAR  code.  A complete  description  of  this  code 
and  its  subroutines  is  contained  in  the  TACWAR  documenta- 
tion1-7** . 

This  study  has  investigated  the  accuracy  of  (1) 
the  nuclear  weapon  effects  information  contained  in  the 
damage  assessment  subroutines,  (2)  selected  algorithms  to 
calculate  target  coverage,  and  (3)  logic  in  the  subroutines. 
The  effort  has  rot  addressed  escalation  criteria,  priority 
criteria  in  nuclear  targeting,  or  the  logic  of  weapon 
assignment  algorithms. 


* TACWAR  was  developed  by  the  Institute  for  Defense 
Analyses . 

**  References  are  listed  immediately  following  Section  II. 
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The  following  subroutines  were  scrutinized  in  this 


'•  NUC6  . The  master  program  for  nuclei!  damage 
assessment. 

• DAMEVL.  The  main  subroutine  which  performs 
damaqe  assessment  calculations  with  the  aid 
of  nine  called  subroutines. 

• DOSLIM.  Sets  the  limiting  nuclear  radiation 
doses  to  do f ini'  "pools"  into  which  personnel 
are  placed,  depending  on  the  dose  they 
receive . 

• FN . Calculates  the  range  to  given  initial 
nuclear  radiation  doses,  the  boundaries  of 
the  pools  established  by  DOSLIM. 

• PREFN’.  A utility  subroutine  which  acts  as 
an  interface  between  FN  and  QK1NR,  and  calls 
QKINR  to  calculate  the  initial  radiation  dose 
at  each  of  23  slant  ranges. 

• QKINR.  Calculates  the  dose  due  to  neutrons, 
air-secondary  gamma  rays  and  f iss ion-product 
gamma  rays  at  each  of  the  ranges  specified 
by  PREFN. 

• WRAD.  Calculates  the  weapon  radius  against 
personnel  in  specified  shelter  postures, 
exposed  to  nuclear  air  bursts  or  surface 
bursts  of  specified  yields. 

• WRADVN.  Implements  the  Physical  Vulnerability 
System3  to  calculate  the  weapon  radius  against 
material  targets  given  the  yield  and  burst 
height  of  the  weapon  and  the  vulnerability 
number  of  the  targec. 

• OFFCOV.  Calculates  the  expected  coverage  of 
a circular  target  of  uniform  value  by  a wea- 
pon delivered  at  an  offset  aimpoint. 

• SIMCN.  Calculates  the  cumulative  circular 
normal  distribution  function. 


10 


r 


The  following  subroutines  were  scrutinized  in  this 

project : 

• NUC6 . The  master  program  for  nuclear  damage 
assessment. 


• DAMEVL.  The  main  subroutine  which  performs 
damage  assessment  calculations  with  the  aid 
of  nine  called  subroutines. 


i- 

F 

F 

| 

i 

i 

i 

i 

i 


• DOSLIM.  Sets  the  limiting  nuclear  radiation 
doses  to  define  "pools"  into  which  personnel 
are  placed,  depending  on  the  dose  they 
receive . 

• FN.  Calculates  the  range  to  given  initial 
nuclear  radiation  doses,  the  boundaries  of 
the  pools  established  by  DOSLIM. 

• PREFN.  A utility  subroutine  which  acts  as 
an  interface  between  FN  and  QKINR,  and  calls 
QK1NR  to  calculate  the  initial  radiation  dose 
at  each  of  23  slant  ranges. 

• QKINR.  Calculates  the  dose  due  to  neutrons, 
air-secondary  gamma  rays  and  fission-product 
gamma  rays  at  each  of  the  ranges  specified 
by  PREFN. 

• WRAD.  Calculates  the  weapon  radius  against 
personnel  in  specified  shelter  postures, 
exposed  to  nuclear  air  bursts  or  surface 
bursts  of  specified  yields. 


• WRADVN.  Implements  the  Physical  Vulnerability 
System8  to  calculate  the  weapon  radius  against 
material  targets  given  the  yield  and  burst 
height  of  the  weapon  and  the  vulnerability 
number  of  the  target. 


• OFFCOV.  Calculates  the  expected  coverage  of 
a circular  target  of  uniform  value  by  a wea- 
pon delivered  at  an  offset  aimpoint. 

• SIMCN.  Calculates  the  cumulative  circular 
normal  distribution  function. 
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SIRCOV.  Calculates  the  expected  coverage  of  a 
circular  Gaussian-distributed  target  by  a weapon 
delivered  at  an  offset  aimpoint. 


• CIRCOV.  Calculates  a Gaussian  function  integrated 
over  an  offset  circle.  It  is  eauivalent  to 
either  the  probability  of  an  offset  aimed  wea- 
pon impacting  a circular  target,  or  that  an  off- 
set aimed  weatxjn  having  a specified  radius  will 
hit  a point  target. 

Detailed  descriptions  of  the  results  of  investiga- 
tion of  these  subroutines  are  presented  in  sections  2-1  to 
2-10  of  this  report.  A summary  of  those  results  is  presented 
below* . 

1-1  DAMEVL 

The  portion  of  DAMEVL  which  assesses  damage  to  per- 
sonnel contains  a call  to  WRADVN  with  a vulnerability  number 
(VN)  for  personnel.  No  source  for  such  a VN  could  be  iden- 
tified. An  alternative  method  using  subroutine  WRAD  would 
require  classified  input  data  and  a rather  subtle  program 
change  to  be  operative.  This  alternative  method  should  be 
added  to  the  program. 

In  addition,  weapon  radii  in  radiation-dominated 
portions  of  WRAD,  which  are  based  on  a 450-rad  lethality 
criterion,  when  compared  with  the  results  of  FN , are  always 
selected  over  those  for  higher  radiation  "states."  As  a 
result,  no  targeted  personnel  would  ever  undergo  transitions 
to  the  higher  states. 

Finally,  DAMEVL  does  not  assess  damage  to  military 
personnel  at  surface  to  surface  missile  sites  or  at  supply 
depots . 

1-2  PREFN 

Opportunity  for  user  input  of  weapon  type,  fission 
fraction  and  air  density  has  not  been  provided.  These  have 


* Note  that  not  all  subroutines  listed  in  Section  II  are 
mentioned  in  this  summary. 
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been  fixed  at  fission,  one-half  and  1.1  x 10"3  gm-cm“3, 

2 

respectively.  Height  of  burst  can  be  either  174  ft/kt3  or 
zero.  Aside  from  contributing  to  unrealistic  input  data, 
this  inflexibility  could  result  in  errors  of  20  to  50  percent 
in  dose  at  a given  slant  range.  An  error  of  a factor  of  100 
exists  in  the  height  of  burst  argument  for  calculating  the 
slant  range,  which  will  cause  significant  errors  in  some 
ground  ranges . 

1-3  QKINR 

Significant  divergences  occur  between  TACWAR  cal- 
culations and  other  newly  calculated  radiation  levels  at 
given  slant  ranges.  A neutron  multicollision  term  (x  2) 
has  been  introduced  into  the  calculations  but  apparently 
not  into  the  shielding  or  biological  response  data.  This 
point  is  the  subject  of  further  investigation.  The  weapon 
types,  "fission,"  "intermediate"  and  thermonuclear"  are  not 
sufficiently  broad  for  force  posture  studies  requiring  the 
capability  to  assess  the  possible  impact  of  enhanced  radi- 
ation weapons. 

1-4  WRAD 

Only  nuclear  radiation  and  air  blast  effects  are 
now  included  in  this  subroutine.  Thermal  radiation  is  the 
dominant  lethal  effect  on  unwarned,  exposed  personnel,  over 
a significant  range  of  yields  from  nuclear  surf  . -e  bursts. 
This  effect  could  be  incorporated  in  the  subroutine.  Yield 
interval  breakpoints,  coefficients  and  exponents  have  not 
been  verified. 
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This  subroutine  now  permits  use  of  only  two  scaled 
heights  of  burst,  zero  and  174  feet.  Consequently,  weapon 
radii  against  soft  targets  (i.e.,  QVN  <_  10,  PVN  < 15)  are 
reduced  up  to  40  percent  because  optiumum  height  of  burst 
cannot  be  utilized. 

1-6  OFFCOV 

The  algorithm  uses  dimensions  normalized  in  units 
of  the  weapon  radius.  However,  prior  to  this  normalization 
a lower  limit  to  the  calculation  is  established  by  setting 
the  coverage  at  zero  if  the  weapon  radius  (WR)  (unnormal- 
ized— any  units)  is  less  than  0.001.  WR  should  first  be 
normalized  in  units  of  target  radius. 

In  the  target  coverage  interval  between  10  and  90 
percent,  the  accuracy  of  the  algorithm  frequently  deviates 
by  as  much  as  20-30  percent  from  the  true  value.  1 n one  ex" 
treme  case,  a numerically  integrated  value  of  0.100  is  calcu- 
lated by  ....  algorithm  as  0.186,  an  86  percent  error.  Mox, 
extreme  values  of  coverage  contain  larger  (percentage) 
errors . 

1-7  SIMCN 

The  algorithm  used  to  replace  the  exact  integra- 
tion is  sufficient. .y  accurate  in  all  relevant  regimes  of  the 
independent  variable.  However,  another  algorithm,  somewhat 
simpler  to  implement  and  slightly  more  accurate,  is 
suggested  for  use.  Numerous  documentation  errors  have 
been  noted  and  tabulated. 
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1-8 


SIRCOV 


This  algorithm  also  normalizes  the  variables  in 
units  of  (adjusted)  weapon  radius.  However,  prior  to  nor- 
malization, the  coverage  is  set  at  zero  for  values  of  the 
weapon  radius  (in  any  units)  less  than  0.1.  A method  of 
normalizing  this  variable,  making  it  dimensionless,  is 
suggested.  Several  documentation  errors  have  been  noted 
and  tabulated. 

1-9  CIRCOV 

The  selected  algorithm  is  accurately  implemented 
in  the  Fortran,  but  has  discontinuities  at  two  points  in 
variable-space,  and  other  inaccuracies  which  produce  10- 
percent  errors  in  predicted  coverage  when  the  coverage  is 
in  the  range  5 percent  to  95  percent.  Outside  this  region 
the  errors  rapidly  increase  to,  for  example,  a factor  of 
two  at  1 percent  coverage. 

1-10  RECOMMENDATIONS 

Based  upon  analysis  of  the  TACWAR  nuclear 
effects  data  and  subroutines,  the  following  changes  or  addi- 
tions to  the  subroutines  in  NUC6  are  recommended: 

a.  Provide  an  algorithm  for  calculating  dose 
from  enhanced  radiation  weapons. 

b.  Allow  only  one  option,  that  employing  sub- 
routine WRAD,  to  calculate  weapon  radius  for  personnel 
damage,  and  set  WRAD  = 0 over  yield  ranges  where  it  is 
radiation-dominated . 

c.  Update  the  dose  calculation  from  initial 
nuclear  radiation  to  conform  to  more  recent,  more  accurate 
calculations . 
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d.  Include  a thermal  option  for  unwarned  person- 
nel exposed  to  a surface  burst. 

e.  Include  personnel  damage  assessment  in 
surface  to  surface  missile  sites  and  supply  depots. 

f.  Implement  various  normalizations  (OFFCOV, 
SIRCOV) , an  algorithm  change  (SIMCN) , and  correct  the 
numerical  error  (SIRCOV)  in  the  several  circular  coverage 
subroutines . 


SECTION  II 


DETAILED  DISCUSSION  OF  SUBROUTINES 


2-1  DAMEVL 

This  is  a multipurpose  subroutine  which  has  among 
its  functions  the  evaluation  of  damage  to  targets  includ- 
ing military  and  civilian  personnel.  As  presently  coded, 
DAMEVL  calculates  personnel  damage  in  the  following 
situations : 


Situation 


Military 

Personnel 


Civilian 

Personnel 


Battlefield  x x 
Airbases  x x 
SSM  Sites*  x 
Supply  Depots*  x 


The  assessment  should  be  extended  to  include  military 
personnel  at  supply  depots  and  at  surface  to  surface 
missile  sites.  Conversations  with  field-experienced  army 
oificers  as  well  as  with  members  of  the  analytical  com- 
munity have  confirmed  our  belief  that  trained  military 
personnel  are  important  to  successful  operation  of  both 
of  these  facilities. 

Radii  of  effects  against  personnel  from  two  dif- 
ferent calculations  are  compared:  (1)  calculations  using 
Part  III  of  reference  8,  and  (2)  subroutine  FN.  These  two 
calculations  are  referred  t^  as  methods  (1)  and  (2)  below. 


* Comparison  of  WRAD  with  FN  is  not  carried  out  in  the 
current  version  of  TACWAR  in  assessing  civilian  per- 
sonnel damage  at  SSM  sites  and  supply  depots. 
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Method  (1)  involves  calculation  of  weapon  radius 
as  in  reference  8 by  use  of  a personnel  vulnerability  num- 
ber (VN)  . No  VN  is  provided,  nor  is  any  acceptable  VN  known 
to  exist.  There  is  another  option,  which  can  be  activated 
by  a change  in  coding.  This  option  relies  on  power-law  fits 
to  the  curves  of  weapon  radius  against  personnel  in  Part  III 
of  reference  8.  For  a given  shelter  posture,  weapon  radius 
against  personnel  depends  on  yield  and  height  of  burst. 

At  low  yields  nuclear  radiation  is  the  dominant  kill  mecha- 
nism and  at  higher  yields  air  blast  becomes  dominant.  Sub- 
routine WRAD  implements  an  algorithm  which  approximates 
the  weapon  radius  vs.  yield  curves  and  calculates  weapon 
radii  against  personnel  in  given  shelter  postures  subject- 
ed to  weapons  of  specified  yields.  This  will  discussed 
in  more  detail  in  Section  2-6. 

Method  (2)  uses  the  initial  radiat  . algorithms 
to  determine  radii  to  the  limiting  dose  levels  of  higher 
radiation  "states"  or  "pools"  (see  discussion  under  sub- 
routine DOSLIM) . Subroutines  FN,  PREFN  and  QKINR  are 
exercised  to  determine  the  radii  from  surface  ground  zero 
to  given  initial  radiation  doses  at  the  boundaries  of  the 
various  "states." 

Weapon  radii  calculated  by  methods  (1)  and  (2) 
are  compared.  The  larger  radius  is  selected  for  calcula- 
tion of  target  coverage.  This  means  that  transitions  into 
higher  radiation  states  can  never  occur  because  for  low 
yield  weapons  the  radius  to  450  rad  will  be  larger  than 
that  for  transitions  to  higher  radiation  states,  and  for 
higher  yield  weapons  blast  or  thermal  radii  will  dominate. 

It  therefore  seems  appropriate  to  set  weapon  radius  to 
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zero  in  the  yield  intervals  of  WRAD  when  the  weapon  radius 
is  dominated  by  nuclear  radiation. 

2-2  DOSLIM 

| This  subroutine  establishes  upper  and  lower 

! boundaries  of  radiation  states.  For  example,  in  the  TACWAR 

! 

i print-outs  examined  in  this  effort  the  following  dose 

i levels  were  selected  as  defining  radiation  states. 

I 

j State  12  345 

j Dose,  rad  0 50  450  3000  8000  >20000 

s 

1 TACWAR  logic  allows  for  statistical  fluctuations  by  placing 

I 

' all  personnel  in  the  center  of  the  dose  range  defining  each 

pool  or  state.  Thus,  exposure  to  25  rad  will  result  in 
transition  from  state  1 to  state  2,  200  rad  w.  il  cause  a 
transition  from  state  2 to  state  3,  etc. 

2-3  FN 

This  subroutine  scans  the  array  >c  doses  calcu- 
lated by  QKINR  as  directed  by  PREFN,  and  findr  the  hori- 
zontal range  at  which  the  dose  is  equal  to  the  input  value, 
which  is  one  of  the  boundaries  of  the  radiation  states 
furnished  by  DOSLIM,  Oi.  t ir  the  case  of  civilian  personnel, 
the  casualty  or  fatality  dose  (250  or  450  rad,  respective- 
ly) . Neutron  and  gamma  transmission  factors  are  first 
applied,  according  to  protection  category  date  for  the 
personnel  being  considered. 

FN  is  called  only  for  military  and  civilian 
targets  at  airbases  and  on  the  battlefield. 
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PREFN 


I 

I 


I 
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This  subroutine  is  used  as  a controlling  sub- 
routine between  FN  (where  range  to  a given  dose  level  is  cal- 
culated) and  QKINR  (where  the  dose  level  at  a given  range  is 
computed) . PREFN  receives  input  parameters  from  DAMEVL  and 
calls  QKINR  at  each  of  23  horizontal  ranges  and  input  burst 
heights  for  every  weapon. 

PREFN  is  used  to  provide  to  QKINR  certain  input  data 
which  set  the  parameters  for  dose  calculations.  These  para- 
meters appear  to  be  frozen  at  single  values,  without  present- 
ing any  options  to  the  user  for  inserting  weapon-specific  or 
situation-specific  values.  Some  inputs  are  in  fact  rather 
unrealistic.  For  example,  only  fission  weapons  and  50-percent 
fission  fractions  are  presently  used.  Only  two  heights  of 
burst,  zero  or  174  ft/kt ' , are  allowed.  Atmospheric  density 
is  set  at  1.1  x 10' 3 gm/cm3.  We  recommend  that  weapon  para- 
meters be  specified  as  part  of  the  input  data  base  and  that 
user  options  be  added  for  air  density  as  well. 

PREFN  also  furnishes  input  on  the  neutron  multi- 
collision factor.  This  factor  is  presently  set  at  2.  The 
effect  is  that  all  free  field  neutron  doses  are  doubled  for 
calculations  of  biological  response.  It  appears  that  the 
factor  is  inserted  correctly  in  defining  the  in-body  radia- 
tion environment,  but  that  it  must  also  be  applied  to  bio- 
logical response  criteria  which  are,  so  far  as  can  be,  based 
on  free-field  doses. 
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One  of  the  functions  of  PREFN  is  to  calculate  slant 
range  to  the  target  and  this  appears  to  contain  an  error.  As 
described  in  the  documentation  IHOB  is  either  1 (airburst)  or 
0 (surface  burst)*  and  if  1,  the  height  of  burst  is  set  at 
1.74,  in  units  of  one  hundred  feet.  Thus  when  calculating 
slant  range  it  appears  that  a factor  of  100  should  multiply 
the  term  HBR  in  order  to  produce  a correct  value  of  slant 
range  to  the  target. 

2-5  QKINR 


This  subroutine,  which  calculates  the  initial 
radiation  environment  from  nuclear  air  and  surface  bursts, 

, 9 

was  verified  against  the  source  material  on  which  it  was 
based,  and  was  also  compared  with  results  of  more  recent 
work 1 c in  the  field  of  radiation  transport  models. 

Calculation  of  dose  from  neutrons  and  air-secondary 
gamma  rays  is  a relatively  straightforward  process  involving 
spherical  divergance  and  exponentials  in  slant  range,  normal- 
ized to  fit  available  empirical  data.  A rather  simple  height- 
of-burst  correction  is  also  employed.  Examples  of  expressions 
used  in  the  QKINR  algorithm  are  contained  in  Appendix  A. 

The  algorithms  used  in  QKINR  for  the  fission  product 
contribution  were  developed  as  a fast-running  implementation 
of  the  method  of  French  and  Mooney9.  They  are  fully  described 
in  a working  paper  provided  by  the  Institute  for  Defense 
Analyses  and  reproduced  as  Appendix  B.  A functional  descrip- 
tion of  the  algorithms  is  also  given  in  Appendix  A. 

Two  routes  were  followed  in  verifying  initial  nu- 
clear radiation  environments  of  TA'  ’ \R: 

• Results  of  independent  ("offline")  use  of  the 
simplified  algorithms  were  compared  with  the 
model  on  which  they  are  based,  and  on  more 
recent  calculations. 
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• QKINR  as  received  from  CCTC  was  exercised 

("online")  to  determine  nuclear  environments 
directly . 

Figures  2-1  through  2-13  contain  results  of  off- 
line calculation  of  the  algorithms  described  in  Appendix  B 
and  the  CCTC  listing.  These  calculations  were  programmed  and 
executed  on  a microcomputer.  The  results  have  been  plotted 
as  total  initial  radiation  dose  as  a function  of  ground 
(horizontal)  range  for  two  types  of  warhead  design,  fission 
and  thermonuclear.  Calculations  were  performed  for  surface 
bursts  and  air  bursts  at  yields  of  1*,  10  and  100  KT.  Com- 
parison standards  were  selected  as  follows.  At  all  three 
yields,  data  from  reference  9 labeled  "French  and  Mooney" 
on  the  figures  are  shown.  For  10  and  100  KT,  more  recent 
calculations10  were  obtained  and  added  to  the  comparison 
base,  labeled  "Gritzner  et  al."  In  addition,  the  10-KT 
thermonuclear  air  burst  was  singled  out  for  more  detailed 
checks  of  the  initial  radiation  components,  i.e.,  neutrons, 
air-secondary  gamma  rays,  and  fission-product  gamma  rays, 
against  available  references.  Data  were  also  obtained,  for 
general  interest,  from  reference  11. 

The  degree  to  which  TACWAR  calculations  agree  with 
comparison  standards  may  be  observed  by  examining  Figures 
2-1  through  2-13.  In  general,  TACWAR  calculations  of  total 
dose  are  in  fairly  good  agreement  with  the  data  of  French 
and  Mooney  on  which  they  were  based.  This  seems  to  be  due 
to  balancing  errors  of  the  dose  components.  TACWAR  is  high 
by  a factor  of  two  in  neutron  dose,  but  low  with  respect  to 
the  fission-product  gamma  dose.  The  total  dose  has  been 

* In  the  interest  of  realism,  only  the  fission  case  was 
run  at  1 KT. 
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Figure  2-1 
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Figure  2-2  Initial  Radiation  Dose  vs  Ground  Range 
for  a 1-KT  Fission  Surface  Burst. 
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Figure  2-4  Initial  Radiation  Dose  vs  Ground  Range 
for  a 10-KT  Fission  Air  Burst  (HOB  = 
114m) . 
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Figure  2-4 
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resolved  into  its  component  parts  (i.e.,  neutrons,  air- 
secondary gammas  and  fission-product  gammas)  for  illustrative 
purposes  in  Figures  2-6  to  2-8. 

Online  calculations  of  these  cases  were  performed 
using  QKINR  on  the  S3  Univac  1108,  and  were  found  to  be  iden- 
tical with  the  offline  calculations. 

None  of  the  existing  weapon  types  in  QKINR  addresses 
the  enhanced  radiation  (ER)  warhead.  The  neutron  spectrum  of 
the  ER  is  sufficiently  different  from  that  of  fission  or 
thermonuclear  weapons  that  a separate  weapon  type  should  be 
included  in  QKINR.  Use  of  ER  warheads  is  being  studied  at 
many  levels  in  the  DOD,  and  it  would  appear  to  be  prudent  to 
include  an  ER  weapon  in  the  calculations  to  allow  inclusion 
of  the  ER  within  the  scope  of  TACWAR  studies. 

2-6  WRAD 

This  subroutine  calculates  the  weapon  radius  due  to 
blast  (or  radiation)  against  personnel  located  in  one  of  eight 
protection  categories.  Input  parameters  consist  of  coeffi- 
cients and  exponents  C and  A in  the  expression 


WRP  = 


3.281 


Where  WRP  is  the  weapon  radius  in  meters  and  W the  yield 
input.  C (in  feet)  and  A are  directly  obtainable  for  various 
protection  categories  as  defined  in  reference  8.  They  are 
not  currently  available  in  unclassified  form,  hence  the  data 
array  for  these  inputs  is  set  to  all  zeros.  The  curves  of 
weapon  radius  vs.  yield  in  reference  8 are  divided  into  por- 
tions of  nuclear  radiation-dominance  and  blast-dominance,  each 
approximated  by  a straight  line  segment,  with  changes  in  slope 
occuring  at  yields  where  radiation  becomes  dominant  over  blast. 
The  yield  at  which  this  occurs  is  called  the  "breakpoint." 

Thus  WRAD  uses  an  algorithm  which  specifies  the  broken  curve 
illustrated  in  Figure  2-14  by  the  constants  C and  A for  each 
region  and  the  yield  breakpoint  separating  the  regions  of 
validity.  ,c 


Lgure  2-14.  Sample  curve  of  weapon  radius  vs  yield,  for  personnel 


There  is  no  allowance  for  thermal  effects  upon 
unwarned  personnel  in  the  open.  Comparison  of  weapon  radii 
due  to  thermal  effects  of  both  "fatal"  and  "incapacitating" 
intensity  shows  that  thermal  fatality  radii  dominate  blast 


fatality  radii  for  nuclear  surface  bursts  at  the  higher  end 
of  the  abscissa.  This  can  be  accounted  for,  if  desired,  by 
adding  or  adjusting  the  exponents  and  coefficients  for  the 
protection  category  for  exposed  personnel. 

Preliminary  investigation  indicates  that  one  break- 
point with  two  coefficients  and  exponents  can  adequately  fit 
the  curves  in  reference  8.  However,  the  actual  values  used 
remain  to  be  fully  verified.  TACWAR  has  the  capability  of 
implementing  a curve-fitting  procedure  using  multiple  break- 
points . 

2-7  WRADVN 

This  subroutine  is  used  to  calculate  the  weapon 
radius  due  to  nuclear  airblast.  Methodology  of  the  Physical 
Vulnerability  System a is  employed.  Input  quantities  are 
the  yield,  height  of  burst  and  vulnerability  number.  Output 
quantities  are  the  weapon  radius  and  damage  "sigma"  as  de- 


fined in  reference  12.  Detailed  investigation  of  the 
Fortran  listing  for  WRADVN  reveals  that  methodology  of  refer- 
ence 8 was  accurately  followed.  Appendix  C contains  the 
Fortran  listing  as  it  appears  in  TACWAR,  with  explanatory 
comments . 


TACWAR  allows  only  two  values  for  the  scaled  height 
of  burst,  zero  and  174  feet.  In  the  case  of  a surface  burst 
the  table  of  weapon  radii,  TABWR,  contains  entries  from 
Table  1-16  of  reference  8 for  P-type  targets  and  from  Table 
1-18  for  Q-type  targets.  For  air  bursts  the  height  of  burst 
(HOB)  inferred  by  interpolating  weapon  radii  was  174  feet, 
in  agreement  with  subroutine  DAMEVL  which  sets  HOB  at  174 
feet  for  all  air  bursts. 


The  vulnerability  number  IVN  is  unpacked  by  the 
usual  methods.  Thus  a VN  of  14Q7  (which  appears  as  1427  in 
TACWAR)  is  resolved  into  its  component  parts:  VN,  the  hard- 
ness; IPQ,  the  target  type  (1  = P,  2 = Q)  and  the  yield- 
adjustment  factor  XK  (or  K in  the  usual  nomenclature  of 
reference  8).  The  "damage  sigma"  (WSIG)  is  set  at  0.2  for 
P-targets  and  at  0.3  for  Q-targets.  Next  the  VN  is  adjusted 
for  yield  in  accord  with  the  methodology  of  reference  8. 

The  resulting  VN  is  used  as  the  entering  argument  for  looking 
up  the  scaled  weapon  radius  in  array  TABWR.  Logarithmic 
interpolation  is  employed  to  correct  for  fractional  values 
of  the  adjusted  VN.  Finally,  the  weapon  radius  is  scaled 
up  (or  down)  to  the  input  value  of  weapon  yield. 

Restricting  all  nuclear  airbursts  to  a fixed 
scaled  altitude  of  174  ft/kt1  has  the  effect  of  reducing  the 
effectiveness  of  weapons  against  soft  targets.  Increases 
of  up  to  40  percent  in  WR  compared  to  the  value  of  WR  at 
174  ft/kt'  can  be  achieved  against  soft  targets  by  increasing 
the  HOB  to  its  optimum  value.  Figure  2-15  shows,  for  P 
type  and  Q type  targets,  respectively,  the  ratio  of  weapon 
radius  at  optimum  HOB  to  that  at  174  ft/kt’  (WRop//WR17  4^  as 
a function  of  the  vulnerability  number  (VN) . Significant 
penalities  in  WR  at  174  ft/kt5  are  encountered  only  for  soft 
targets.  Table  2-1  illustrates  this  for  VNs  represent- 
ing soft  targets.  Thus,  for  P VNs  of  about  15  or 
less,  or  Q VNs  of  20  or  less,  the  use  of  optimum  HOB  becomes 
important  from  the  standpoint  of  increased  WR.  This  situation 
can  be  corrected  by  adding  3 additional  scaled  heights  of 
burst,  i.e.,  six  aaditinal  tables  in  array  TABWR,  and  suit- 
able programming  to  allow  matching  of  weapon  HOB  to  target 
hardness.  Moreover,  direct  consultation  with  members  of  the 
IDA  staff  indicates  that  memory  limitations  do  not  preclude 
a more  detailed  listing  of  WR  as  a function  of  HOB. 
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Figure  2-15  Increase  in  Weapon  Radius  by  Use 
of  Optimum  HOB 
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Table  2-1  Comparison  of  weapon  radii  at  optimum  air  burst 
altitude  with  those  at  the  standard  height  used 
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For  targets  harder  than  15P  and  10Q,  the  relative 
gain  in  weapon  radius  at  the  optimum  HOB,  compared  to  that 
at  174  ft/kt1,  is  diminishing  rapidly.  Accordingly,  a con- 
stant value  of  WR  for  0 <_  SHOB  < 400  and  VN  >10Q  (or  15P) 
would  be  satisfactory. 

There  is  another  aspect  of  the  height-of-burst 
problem  which  should  be  considered.  In  some  wargaming 
studies,  it  may  become  necessary  to  depart  from  unclassified 
representative  yields  of  tactical  and  theater  weapons  and  to 
assess  the  impact  of  the  use  of  real,  operational  weaponry 
for  both  sides.  To  accommodate  such  a need  it  will  be  nec- 
essary to  extend  the  height-of-burst  argument  to  the  entire 

l 

range  from  zero  to  900  ft/kt’,  and  to  add  appropriate  program 
steps.  For  example  a given  yield  and  height  of  burst  might 
result  in  a scaled  height  of  burst  (SHOB)  of  650  ft/kt'.  In 
this  case  linear  interpolation  would  be  required.  For  a 

i 

SHOB  of  400  ft/kt’,  a weapon  radius  between  those  for  174 
and  600  ft/kt1  will  be  calculated  by  interpolation. 

2-8  OFFCOV 

This  subroutine  calculates  the  expected  coverage 
of  a uniform-value  circular  target  by  a cookie-cutter  weapon 
(that  is,  one  having  a distance-damage  sigma  of  zero)  aimed 
with  Gaussian  aiming  arrors  at  an  offset  aiming  point.  Thus, 
the  independent  parameters  are  the  CEP  of  the  weapon,  the 
target  radius,  the  offset  distance,  and  the  weapon  radius. 

Any  self-consistent  set  of  units  may  be  used,  but  in  fact 
the  calculations  are  made  in  units  of  weapon  radius. 
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The  Fortran  listings  as  provided  by  IDA  have  been 
examined  for  consistency  and  logic  flow  (to  confirm  that  all 
areas  of  input  variable  phase-space  have  been  covered) , but 
have  not  been  verified  against  any  original  documentation. 

The  only  possible  problem  identified  is  at  the  beginning  of 
the  subroutine  in  which,  prior  to  normalizing  all  dimensions 
into  units  of  weapon  radius,  the  coverage  is  taken  to  be  zero 
for  values  of  the  weapon  radius  WAN  < .001.  This  will  lead 
to  a gross  error  if  the  units  of  dimensions  are  large,  and 
will  lead  to  unnecessary  calculations  through  the  main  algor- 
ithm if  the  units  are  small,  such  as  meters.  Normalization 
in  units  of  the  target  radius  (TER)  would  be  appropriate: 
if  WAN/TER  < .001,  set  coverage  equal  to  zero. 

A complete  flowchart  complementing  that  shown  in 
figure  82  of  the  CCTC  documentation13  was  developed  and  pro- 
grammed on  a Tl-59 , with  the  exception  of  the  called  sub- 
routine CIRCOV  used  for  normalized  values  of  the  target 
radius  less  than  0.2.  The  results  of  these  calculations 
using  the  OFFCOV  algorithms  are  compared  to  the  results  ob- 
tained by  direct  numerical  integration  (kindly  provided  by 
L.  Schmidt  of  IDA)  in  Appendix  D,  Tables  D-l  to  D-33.  Illus- 
trative samples  of  these  data  are  shown  in  Figures  2-16 
through  2-21  for  values  of  the  normalized  target  radius  in 
the  range  0.2  <_  TAR  < 2 , in  which  the  error  in  the  coverage 
function  COV  is  plotted  against  the  normalized  aimpoint  off- 
set for  parametric  values  of  the  CEP.  Errors  do  not  exceed 
30  percent  for  values  of  the  independent  variables  which 
lead  to  target  coverages  in  the  range  0.1  <_  COV  0.9. 

Users  of  the  TACWAR  model  should  be  aware,  however,  that  if 
extreme  targeting  regimes  are  being  used  (such  as  preclusion 
oriented  targeting  in  which  a collateral  damage  "target" 
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Figure  2-20.  OFFCOV  - Error  in  COV  for  TAR 


coverage  <0.1  might  be  required,  or  in  which  there  are  very  high 
coverage  requirements  of  >0.9)  that  the  use  of  OFFCOV  can  lead 
to  large  errors.  The  largest  absolute  error  identified  is  for 
TAR  = 0.2,  CEP  = .04,  OFFSET  = 0.75,  for  which  the  OFFCOV 
algorithm  calculates  COV  = 0.848,  whereas  the  numerical 
integration  gives  the  value  of  COV  = 0.688. 

The  documentation  2 ' 5 ' 6 was  not  specifically 
reviewed  for  errors,  but,  where  errors  were  noted,  they  have 
been  included  here.  Since  these  errors  have  in  general  not 
been  carried  over  into  the  Fortran  statement,  they  will  not 
lead  to  any  error  through  use  of  the  TACWAR  model  itself, 
but  could  lead  to  some  confusion. 


1.  Reference  2,  page  C-4 , the  constant  1.1744 
was  used  rather  than  1.1774. 


2.  Reference  6,  page  490,  the  lower  left  box  of 


flowchart  contains  the 

expression : 

COV  = PNO  * exp 

/ OFFSET]2! 

- \ 2o  / J 

COV  = PNO  * exp 

(OFFSET)2! 

2 (o  ' ) 2 J 

which  should  be 


2-9 


SIMCN 


This  subroutine  calculates  the  cumulative  normal 
distribution  function: 


By  defining  x = j ARG/^2~  | , the  integral 
form  of  the  error  function: 


«] 

can  be  put  in  the 
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ARG  > 0 
ARG  < 0 

*"  O 

An  approximation11*  is  used  to  evaluate  erf  x (defined  as  the 
variable  CUP  (x)  in  SIMCN) : 

erf  x=l  - (1  + ax  + ax2  + ax3+a  x4)-4, 

12  3 >♦ 

where  a = 0.278393  a = 0.230389  a = 0.000972  a = 0.078108 

12  3 4 

and  having  a stated  error  e,  such  that|c|£  5 X 10  4 . Thus, 

P (ARG)  = H [1  + CUP  (x)  ] , (1) 

defined  as  the  variable  COV.  The  Fortran  listing  for  SIMCN  as 
supplied  by  IDA  correctly  implements  this  algorithm. 

An  alternative  approximation15  can  be  used  for  P(x), 
without  converting  to  the  error  function  first: 

P(x)  = 1 - S(1  +cx+cx2+cx3+c  x4)-4,  (2) 

12  3 4 

where  c = 0.196854,  c = 0.115194,  c = 0.000344,  c = 0.019527, 

12  3 4 

and  having  a stated  error  e,  such  that|cj<^  2.5  X 10  4 . Table  2-2 
compares  the  tabulated  data  1 6 with  the  results  obtained  from 

algorithms  (1)  and  (2)  over  a range  of  values  of  ARG.  Nega- 
tive values  of  ARG  are  not  included  since  they  will  have 

the  same  accuracies.  Note  that  algorithms  (1)  and  (2)  are 
essentially  identical  for  the  trial  values  tabulated,  dif- 
fering by  about  one  in  the  seventh  decimal  place.  Neither 
of  these  algorithms  differs  from  the  exact  tabulated  data  in 
the  NBS  handbook  by  more  than  about  two  in  the  fourth  decimal 
place.  Thus,  either  is  totally  adequate  for  the  purposes  of 
TACWAR . 


P (ARG)  = h [l  + erf  x]  if 
= >5  [l  - erf  x]  if 


x 


in  which  erf  x 


.*/ 


-t* 


dt 
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Numerous  errors  exist  in  the  documentation  we 
have  available.  The  CCTC  description 17  of  SIMCN  is  in 
error  in  its  last  two  expressions  on  page  181.  They 
should  read: 

COV  = P(ARG)  = (Jj)*(l  + CUP  (ARG/ 
if  ARG  > 0 

COV  = P (ARG ) = ( 4 ) * ( 1 - CUP  (ARG/  ^2)) 
if  ARG  < 0 

The  IDA  documentation18  has  errors  in  five  expressions 
under  the  algorithm  implementation.  The  correct  expressions 
are : 


-fc2/2  dt 

<Mx)  = h[l  - 4>(x)]'  if  ARG  < 0 
ip  (x)  = *5  [ 1 + 4>  (x)]  if  ARG  > 0 

The  Fortran  listing  itself  contains  none  of  these 
errors,  so  they  are  not  carried  forward  in  the  model. 


A = 0.078108 


x 


— oc 
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2-10 


SIRCOV 


This  subroutine  calculates  the  expected  coverage  of 
a Gaussian  circular  target  with  offset  aimpoint  by  a weapon 
having  Gaussian  aiming  errors  and  a cumulative  log-normal 
distance-damage  function.  Thus,  the  variables  used  are  the 
aimpoint  offset,  the  weapon  radius,  the  CEP,  the  target 
radius  (containing  95  percent  of  the  target  area) , and  the 
distance-damage  sigma  (in  units  of  the  weapon  radius) . The 
algorithm  implemented  is  that  developed  by  Mason19. 

The  Fortran  listings  as  provided  by  IDA  have  been 
compared  to  the  referenced  document  and  reviewed  for  obvious 
errors.  At  tne  beginning  of  the  subroutine,  prior  to  nor- 
malizing all  dimensions  into  units  of  (adjusted)  weapon 
radius,  the  coverage  is  set  to  zero  for  values  of  weapon 
radius  WPN  < 0.1.  Since  it  is  stated  that  the  variables  may 
be  expressed  in  any  consistent  units,  a large  error  may 
result  if  the  units  are  large,  while  unnecessary  calcula- 
tions will  be  made  through  the  main  algorithm  if  the  units 
are  small.  Use  of  the  variable  R (R  = WPN/SIG) , defined 
later  in  the  statement,  might  be  appropriate. 

There  is  also  a Fortran  statement  error  in  defining 
the  adjusted  CEP,  taken  to  be: 

CEPP  = SQRT  (C*C+0 . 21 3*T*T) . 

The  coefficient  of  T2  is,  in  fact,  the  ratio  In2/ln20, 
which  is: 

In2/ln20  = 0.231378. 

This  digit  inversion  also  appears  in  reference  19.  it  will 
cause  an  error  of  n,  4 percent  at  most  in  the  adjusted  radius 
(for  those  cases  in  which  the  CEP  is  small  compared  to  the 
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target  radius) , and  is  used  only  in  those  parts  of  the  sub- 
routine in  which  the  weapon  distance-damage  sigma  is  either 
0.4  or  0.5.  The  minor  impact  which  results  over  a range  of 
variables  will  be  discussed  later. 

An  exact  (to  the  fourth  decimal  place)  numerical 
integration  of  this  coverage  function  is  available"1’. 
Normalization  is  in  units  of  the  weapon  radius  divided  by 
the  adjusted  CEP  (WR/CEP  ) , and  the  offset  distance  divided 
by  the  weapon  radius  (D/WR) . Tabular  data  are  presented  for 
distance-damage  sigma  of  0.1  to  0.5  in  steps  of  0.1. 

Selected  data  from  this  report  are  shown  in  the  first  row 
of  each  double  row  in  Tables  2-4  through  2-8. 

A detailed  flow  chart  was  made  of  the  SIRCOV  algo- 
rithm and  programmed  on  a TI-59.  The  results  sampled  over 
those  regions  of  variable  space  which  produce  a target  cover- 
age generally  in  the  range  0.1  - 0.9  are  shown  in  the  second 
row  of  each  set  in  Tables  2-3  through  2-7.  It  can  be  seen 
that  errors  in  the  coverage  are  generally  less  than  15  per- 
cent, and  become  this  large  only  when  the  coverage  drops  to 
the  0.1  extreme.  As  an  alternative  to  the  tabulation  of 
the  fractional  error  of  the  fractional  coverage,  Figures 
2-22  through  2-2t  show  curves  of  the  error  in  edveraqe  for 
the  parameters  WR/CEP  , D/WR,  and  sigma  They  generally 
fall  in  the  range  +0.07  > A COV  N -0.14,  with  the  algorithms 
for  SIGT  =0.3  and  0.5  being  the  least  accurate. 

Tables  2-8  and  2-9,  using  the  same  normalization  of 

variables  as  above,  show  the  error  caused  in  the  fractional 
coverage  arising  from  using  the  incorrect  constant  0.213  in 
calculating  the  adjusted  CEP,  as  discussed  above.  The  first 
row  of  each  group  uses  the  correct  constant  0.231  in  the 


54 


\ 

* 


Table  2-'i.  OFFSET  CIRCLE  PROBABILITIES 
Distance  Damage  Sigma  (SIGT)  =0.1 
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Table  2-4.  OFFSET  CIRCLE  PROBABILITIES 
Distance  Damage  Sigma  (SIGT)  = 0.2 
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Table  2-9.  OFFSET  CIRCLE  PROBABILITIES 
Error  Arising  From  Use  of  Constant  0.213. 
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SIRCOV  algorithm,  the  second  the  incorrect  one.  The  largest 
error  found  was  ACOV  = 0.01. 

Although  the  documentation  was  not  specifically 
examined  for  errors,  the  following  were  noted.  They  have 
not  in  general  been  carried  over  into  the  Fortran  statements. 


read : 


The  correct  expressions  in  reference  2 should 


Page  C-8,  last  line:  CEP'  = [c  + 0.231  (T)  2]  ** 

Page  C-8,  last  line:  If  SIGT  < 0.301... 

Page  C-9,  line  11:  P(?.0'  = CUMN 

r,  WPN  ( 1- SIGT 2 ) / 

In  5 / . - etc  . 

! CEP1  / 


Page  C-9,  last  lint:  ...and  r > 1.5, 


5.  Page  C-10,  line  3:  PD  = P(O)  exp  '-|^P(0)r2/R 

On  page  492,  reference  6,  the  flow  chart  reflects 
the  unit  normalization  problem  discussed  above  in  connection 
with  the  Fortran  statement  in  which  WPN  ^>0.1  is  used  rather 
than  the  suggested  R >_  0.1. 
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CIRCOV 


This  subroutine  calculates  the  circular  coverage 
function,  which  is  a Gaussian  function  integrated  over  an 
offset  circle.  It  is  equivalent  to  either: 

a.  Probability  that  an  impact  point  will  fall 
within  a circle  of  radius  R aimed  with  an  offset  r if  the 
impact  has  a Gaussian  probability  of  unit  standard  deviation. 

b.  Probability  that  a point  target  will  be  covered 
by  a weapon  having  an  effects  radius  R (and  a distance- 
damage  o = 0,  i.e.,  a "cookie  cutter"  weapon)  when  aimed 
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I 


I 

f 


A 


with  an  offset  r with  an  aiming  error  having  a Gaussian 
probability  of  unit  standard  deviation. 


The  integral  to  be  evaluated  is: 
P(R.r)  = j + r’’/: 


dx  dy, 


in  which  the  integration  is  carried  out  over  the  circle 
(x-a) 2 + (y-b) 2 = R2 


where 


a + b‘  = r " 


In  practice,  when  the  distances  (including  the  CEP 
aiming  errors,  which  are  generally  not  unity)  are  given  in 
any  consistent  set  of  units,  they  are  normalized  by  dividing 
by  (CEP/1 . 1774)  . 

In  TACWAR , this  integral  is  evaluated  by  an  algor- 
rithm  using  an  approximation  in  the  reference  14  equations 
26.3.25,  26.3.26  and  26.3.27.  The  Fortran  statement  has  been 
carefully  checked  and  found  to  correctly  implement  the  algo- 
rithm. It  was  programmed  on  a TI-59  to  generate  the  data 
shown  in  Figures  2-27  and  2-28,  in  which  discontinuities  at 
R = 1,  5 car.  be  noted.  To  check  the  accuracy  of  the  algo- 
rithm, a numerical  integration  of  P(R,r)  was  made  for  a 
series  of  parametric  values  of  the  independent  variables.* 
These  (having  an  accuracy  of  several  in  the  fifth  digit)  are 
shown  in  Table  2-10,  along  with  the  corresponding  values 
from  CIRCOV,  and  are  also  shown  as  points  on  Figures  2-27 
and  2-28.  it  can  be  seen  that  for  probabilities  in  the 
range  0.05  -0.95  the  CIRCOV  algorithm  is  generally  accurate 
within  10  percent.  It  should  be  kept  in  mind,  however,  that 


i 


i 

I 


* These  data  cover  a broad  spread  of  parameter  space.  In 
some  regions  they  have  been  confirm  : by  comparison  to 
more  detailed  tabulations  such  as  inference  23. 


* 
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for  extreme  probabilities  these  approximation  errors  increase 
to  about  a factor  of  two  at  P = IQ-2  and  a factor  of  five  at 
P = 10“5. 

There  are  several  minor  documentation  errors,  not 
carried  over  into  the  Fortran  and  hence  having  no  impact  on 
the  proper  operation  of  the  subroutine.  References  2 and 
5 fail  to  note  that  for  R > 5 and  RLS  < 1,  then  ARG  (for  the 
called  subroutine  SIMCN)  is  set  equal  to  R,  although  this 
branch  is  correctly  noted  on  Figure  85,  reference  6.  Other 
small  errors  are: 


1.  Figure  85,  reference  6,  the  expression 


is  incorrectly  written. 


and 


2.  Page  185,  reference  5 has  the  inequalities 
R > 1 

1 < R <_  5 incorrectly  written. 
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APPENDIX  A 


FUNCTIONAL  DESCRIPTION  OF  ALGORITHMS 
EMPLOYED  BY  SUBROUTINE  QKINR 
FOR  CALCULATING  INITIAL  RADIATION  DOSES 


The  following  algorithms  are  used  in  the  QKINR 
subroutine  to  calculate  components  of  total  dose  from 
initial  nuclear  radiation.  This  information  has  been 
taken  from  program  listings  and  from  the  IDA  paper  re- 
printed in  Appendix  B. 


Neutron  Dose 


D = 

n 


fCDB  + ADB  • R 2 + BDB  • R 1 
n J 


where  Rq  is  the  slant  range  from  the  detonation  to  the  target, 
TF  is  a height-of-burst  correction  (1.0  for  surface  burst  and 
a fixed  value  for  an  air  burst.  The  value  of  TF  is  dependent 
only  on  the  weapon  type.),  CF  the  neutron  multicollision  fac- 
tor, RBE  the  neutron  relative  biological  effectiveness,  and 
W is  the  total  yield  in  kt.  Parameters  CDB , ADB,  etc.,  and 
TF  are  defined  in  Table  A-l. 

Air-Secondary  Gamma  Dose 


2SG  " 


TF  [CDA  + ADA  • R 2 + BDA 
e ° 


Variable  definitions  are  similar  to  those  given  above. 


Table  A-l.  Values  of  Exponents  and 
Coefficients  Used  in  Calculating  Neutron  and 
Air-Secondary  Gamma  Ray  Doses  in  Subroutine  QKINR 


Parameter 

Fission 

Weapon  Type 
Intermediate 

Thermonuclear 

ADA 

3.67E-8* 

1.00E-7 

1. 28E-7 

ADB 

3.48E-8 

9 . 26E-8 

9 . 36E-8 

BDA 

-4.42E-3 

-5.09E-3 

-5. 15E-3 

BDB 

-3.22E-3 

-3. 81E-3 

-3.70E-3 

CDA 

22.46 

22.91 

21.95 

CDB 

20.67 

20.62 

19.73 

TF 

2.16 

2 . 20 

2.19 

* Read  as 

3.67  x lO-8 

Fission-Product  Gamma  Ray  Dose 

Calculation  of  the  prompt  dose  from  gamma  rays 
emitted  by  mixed  fission  products  is  complicated  by  the 
temporal  and  spatial  dependence  of  the  source  spectrum, 
medium  absorption  properties  and  source  geometry.  The 
following  expressions  were  used  in  the  calculation  of 
fission-product  gamma  ray  dose  from  1 kt,  10  kt  and  100 
kt  surface  burst  and  air  burst  weapons.  They  were  derived 
from  Appendix  B,  which  contains  complete  instructions  and 
algorithms  for  use  with  any  set  of  parameters  of  interest. 

Tpp  q p^c  _ o 5 5 . 2 

1 kt  surface  burst:  D__  = — lO0-0*3  660  R 

FP  p 2 O 

H0B=lm,  Rq  > 362m**  o 


**  If  the  range  decreases  below  a minimum  range  parameter 
(defined  in  Appendix  B)  additional  terms  are  introduced 
in  the  algorithm.  Representative  values  of  minimum 
range  are  362m,  415m,  and  563m  for  1,  10,  and  100  kt, 
respectively.  Ranges  less  than  these  are  generally  not 
of  great  interest  because  of  (a)  extremely  high  dose 
levels  and  (b)  very  high  neutron  dose  relative  to  gamma 
ray  doses. 
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A. 


INTRODUCTION 


A computer  subi'ou t. Inc  l'or  the  calculation  of  the  Initial 
Nuclear  Radiation  Dose  has  been  developed  by  C.  Fisonhauer  and 
L.  Spence  of  the  National  Bureau  of  Standards.  This  subroutine 
computes  the  dose  from  fission  product  gamma  radiation,  secondary 
gamma  radiation,  and  neutrons.  The  calculations  are  based  upon 
a paper  by  French  and  Mooney.1  The  calculations  for  secondary 
gamma  rays  and  neutrons  arc  made  from  equations  of  the  form 

D = a exp  [b  exp  ((cR^)'")  + d exp  (cR  )]  , 

where  R is  the  slant  range  from  weapon  to  monitor  point.  This 
can  he  implemented  for  relatively  rapid  calculation.  The  calcu- 
lations for  the  fission  product,  gamma  dose,  on  the  other  hand, 
are  much  more  lengthy.  They  require  more  complex  expressions 
for  dose  ns  a function  of  time.  Those  expressions  must  be 
integrated  as  a function  of  time  to  obtain  the  total  dose.2  As 
a result  the  computer  routine  which  implements  these  equations 
is  a relatively  slow  calculation.  l'or  analysis  where  the  doses 
must  be  calculated  a largo  number  of  times,  in  particular  in 
damage  assessment  calculations  involving  many  weapons  and  monitor 
points,  such  long  running  calculations  can  add  most  substantially 
to  computer  requirements. 

An  algorithm  for  a more  rapid  calculation  of  the  fission 
product  gamma  doso  is  described  below.  This  algorithm  is  a strictly 
numerical  fit  to  the  results  of  calculations  of  fission  product 
gamma  doses  over  a range  of  parametric  values.  The  fit  is  gener- 
ally within  10  to  20  percent  over  the  range  of  interest.  This 
range  is  for  yields,  W,  ranging  from  0.1  t.o  10,000  KT,  scaled 
heights  of  burst  (the  height  of  burst/W'' ^ ) from  0 to  1,300  feet, 

1 R.  L.  French  and  L.  G.  Mooney,  "Initial  Radiation  Exposure  from  Nuclear 
Weapons,"  Radiation  Research  Associates,  Inc.,  Interim  Report  on  OCD  Contract 
No.  DAHC20-72-C0123 , RRA-T7201,  15  July  1972. 

2Thc  essential  reason  for  the  integration  is  due  t.o  the  hydrodynamic  buoyancy 
of  the  fireball  which  not  only  changes  tlx?  distance  between  source  and  receiver, 
but  due  to  changes  in  air  density  with  height  changes  tlx?  radiation  ifllsorption 
in  a complex  manner. 
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and  slant  ranges  ranging  from  either  the  minimum  slant  range 
at  ground  zero,  or  the  minimum  slant,  range  'where  the  fission 
product  gamma  dose  is  less  than  1/20  the  neutron  dose,  to  a 
maximum  slant  range  wnere  the  fission  product  gamma  dose  is 
about  10R.  The  algorithm  will  be  described  directly  first, 
followed  by  a few  comments  concerning  its  development. 


B.  ALGORITHM  FOR  FISSION  PRODUCT  GAMMA  DOSES 
A maximum  slant  range,  SRM,  is  computed  by 

SRM  = 1 , 8 Q S + 8A8.25-L  + 111  . 2 -’t  1 - L2  , 

where  L = log1Q(W).  If  Ro  > SRM , the  Dose,  1),  is  0;  otherwise 
the  fallowing  procedure  is  followed. 

An  "asymptotic  logarithm  I r " dose,  , i5  computed  as 

follows.  Let 

R„  = 1,188.13  + 8 A . 3 2 b 9 L + 71  . * U0  3L‘  + 12.A08AL3  - U.7729L 
+ 2 . 837  AlL^  . 


D 


asy 


7 


(R  -R  )/(:0Q  . 

o 7. 


If  Hb  > 0.2,  set  Dasy 


.0  + R.  , 

asy  A 


1 /I 

where  HQ  = scaled  height  of  burst  ((ft/(KT)  J)/100), 

RA  = + C*HB  ’ 


with 


b = 1.?A15A*10~2  + 6.0937*10_3L  + 3 . A6bAb • 10~3L2 


+ 1.53 A *10  3L3  - 5.9337  * 10"C,L14 


c = 9.9037*  10"5  9 . 1 A 7 • 1 0 5I.  + l . 963  • lO^L2 

+ 1.836l6*10“\3  - 1. 0b9*  10-/,L^  + 1 . 8] 62 • 10~5L 

Compute  a difference  dose,  Dp,  by  the  following  procedure. 

If  L < 2. A: 

S0  = 362  + 7 A . 3 • L - 55-99-L2  + 3A.59-L3  . 
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K<'r  K n 

o M ’ 


!V  " Do  ( M1<x“Xo)  , 

wln'iv  X = IOO/K  ; 

o 

Do  = -0.01b  - 0.0056Mb  ; 
xQ  = 0.055  - 0.00135Hb  ; 

m = "’o  + mSHB  • 
and  for  L < j 

rn0  - 0-552  + 0.398L  + o.25L?  j 

mS  = °-°5l8  + 0.02985L  + D.01685L?  ; 
and  for  I,  ->  i 

mQ  = -1.71  + 2.78L  ; 

mS  = 0.1690  - 0.0615L  . 


F°r  Ro  < sb»  lGt 

DP  = Do  + m<x-x0)  + c(x-100/SD)2  , 
where  Dq,  xq,  and  m are  computed  as  before  and 


, JO. 22  + [ 0 . 5 7 ( L+ 1 ) ] ^ , L < 1.602  , 
'0.22  + d.mJKL-O.Ol)]'1  , L > 1.6 


602  . 


Now  if  l > 2d,  we  have: 


lf  Ro  i SB 


Dp  = ,n( x-x  ) , 


where  x = 100/R 

o 


m - -7/1.8819  + 95.33871.  - 80.2997L2  * 6.06288L3  ; 
xo  = 0 . 0217  ( '1 . 76-1.)  . 

But  if  Ro  < sB>  compute 

V “ « “ eRo  » when  R0  > 200  , 

D0  a 8R0  + i(200-Ro)  , when  R < 200  , 


’ W ■<:  *h. 


> --V.  - it  ;5«tSi»4fr  :n*l  W ■■ 


15 

iS 


-0  . ) ‘I t (I . ! I i) . o | ; 

)o.ooi«.’3  - o.odo'kui,  , i,  s.  <i  , 

)().  UOlVOfi  - (kOOUlFtfO,-^)  , I.  > i|  ; 
i . • I o"''  - *..(k,8-  L0~  I,  + 'l.Vc’jvlO'  *!/•  . 


Now  lot. 


|V  ■ "b  i.  > 

Do  * + n!i*HB  ' "b  > °-‘J  • 


whore  D„  = O.OjBw  - 0.01!:,  . 


And  let 


Finally,  let. 


and 


1)  = 10 


D = D - D , 
r any  F 


D. 


D " = 


W - 10 


R. 


and 


D = D 


!•’ 


0.5  ’ 

where  F is  the  weapon  fission  fraction, 


C.  COMMENT 

The  underlying  motivation  of  the  above  schema  war.  obtained 
by  observing  that  the  dose  as  a function  of  Ro  asymptotically 
approaches  an  expression  of  the  form 

A*Wexp( -R  ) 


which  would  be  obtained  from  a point  source  with  no  fireball  rise 
arid  constant  absorbing  cross  section.  The  dose  becomes  close  to 
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for* 


till:'  l:'.VIi;pt  Ol  ! I'  ;i!  tin;-.,-  level;-,  e 

'I’lui;-  the  first,  effort  l t.o  obtain  a .linear 


100  t.o  1 , OnoH  . 

It  at  far  ruiij'es 


This  was  done  a 
curves.  The  fit,  was. 
height.  o f burst  was 
ranged  from  about  30 


osumLng  * lie  same  slope  for  all  asymptotic 
first,  made  with  !ln  - 0;  a correction  for 
then  added.  The  height  of  burst:,  correct  lor 
percent  (at  maximum  height  of  burst)  for 


0,1  KT  yields  to  somewhat  over  at  largo  yields. 


hs.tnr,  the  "asymptotic  dos.e,"  the  logarithm  of  the  ratio  of 
asymptotic  to  act.ua  1 dose,  l>p,  was  es.t  i mated.  This  Is  a function 
'■hat  has  high  values,  lor  low  slant,  ramies  and  deerons.es  to  sere 
us  the  two  doses  approach  each  other.  Dp  as  a function  of  i/R  ^ is 
almost  linear  near  the  origin,  followed  by  a segment  which,  for 
most  yields,  was  approximated  by  a parabola  that  is.  tangent  to 
tlie  linear  piece  at.  their  intersection.  The  intersection  occurs, 
where  the  slant  range  has  a value,  bid.,  that  was  determined  by 
inspection  from  granhs  of  the  function.  At  a particular  height 
of  curst. , the  linear  segments  for  all  yields,  below  2'jQ  KT  could 
be  taken,  without  too  much  forcing,  to  have  one  common  Intersection, 
tor  larger  yields  to  have  another  different  common  intersection. 

I ii l s naturally  sopa ra tod  the  calculations,  into  two  ramies  of  yields, 
below  <?r)0  KT  and  above  fhO  KT.  These  intersections  were  height  of 
burst  dependent  for  low  yields,  but  could  be  taken  as  constant  for 
iilflh  yields.  ('Lite  ordinate  of  the  Intersection  is  negative,  which 
results  from  errors  in  est  imat  l nr.  the  asymptotes.  In  effect,  the 
estimation  of  Dp  also  partially  compensates,  for  errors  in  the 
asymptote  estimate,  and  gives  a two-step  correction.) 


1 he  slope  ol  the  linear  sections  was.  represented  by  a linear 
function  of  height  of  burst  for  low  yields,  with  the  coeff  1 0 i ents. 
for  the  linear  function  yield  dependent.  For  the  high  yields,  no 
height  of  burst  sensitivity  was  needed. 

For  small  ranges,  large  values  of  l/»n,  a parabolic  segment 
was  added  to  the  linear  variation  whose  coefficient  was  yield 
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dependent  in  the  low  yield  range.  For  the  large  yield  range , 
thin  procedure  cave  an  inadequate  fit,  so  Tor  low  values  of  II 

o 

an  alternative  procedure  was  used,  namely  estimating  log]0(Dr)  a3  a 
function  of  R0>  A linear  function  was  adequate,  except  for  values 
of  0^  under  200  feet,  where  a parabolic  segment  was  added. 

The  algorithm  used  may  seem  a rather  jerry-built  assemblage 
of  curve  fitting  procedures,  as  in  one  sense  it  is.  The  numerical 
values  were  obtained  either  from  graph  paper  or  simple  least 
squares  polynomial  fits.  The  rationale  for  this  approach  is  that 
a function  of  three  variables  is  to  be  fit,  and  there  is  no  a 
priori  way  of  determining  the  functional  forms  needed  for  efficient 
fitting.  The  variation  of  dose  as  a function  of  slant  range  was, 
in  fact,  well  approximated  as  a ratio  of  two  polynomials.  Unfor- 
tunately the  coefficients  of  these  polynomials  did  not  systematically 
vary  as  a function  of  yield,  or  height  of  burst,  rendering  the 
development  of  an  approximation  valid  for  any  yield  or  height  of 
burst  difficult.  A simultaneous  estimation  technique  with  all 
three  Independent  variables  Included  seemed  required.  Although 
this  was  not  attempted,  it  appeared  likely  that  rather  high  order 
terms  would  be  needed  for  any  adequate  polynomial  approximation. 

Thus  the  method  of  "cut  and  fit"  seemed  more  appropriate. 

The  original  algorithm  and  the  approximation  were  implemented 
On  a Control  Data  6*100  computer,  and  compared  over  a range  of 
yields,  heights  of  burst,  and  slant  ranges.  The  average  time  per 
calculation  of  ail  three  tvnes  of  doses  for  the  original  algo- 
rithm was  0.6*10  seconds,  and  for  the  approximation  0.00176  seconds. 

A display  of  the  accuracy  of  the  approximation  is  presented 
in  Table  1 where  the  minimum  arid  maximum  values  of  the  ratio  of 
fission  product  of  doses  computed  by  the  approximation  to  that,  com- 
puted by  the  numerical  Integration  Is  presented  for  various  yields 
and  scaled  heights  of  burst  over  slant  rang.es  of  interest.  The 
slant  range  of  interest  for  this  table  is  defined  as  any  slant 
range  where  the  fission  product  gamma  dose  is  over  10R,  and  where 
20  times  the  maximum  or  either  the  estimated  or  actual  fission 
product  gamma  dose  is  less  than  the  neutron  dose.  As  can  be  seen, 
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Table  1.  MINIMUM  AND  MAXIMUM  RATIOS  OF  ESTIMATED 
TARGET  DOSES 


Scaled  Height  of  Burst 

(ft/(KT)1/3) 

Yield 
( KT ) 

'0 

100 

180 

400 

750 

1250 

0.1 

0.99 

0.99 

0.99 

1.00 

1.00 

1.00 

0.99 

0.99 

0.99 

1 .00 

1.00 

1.00 

1 

0.98 

0.98 

0.99 

0.99 

0.99 

0.99 

0.99 

1 .00 

1.00 

1.00 

1.00 

1.00 

10 

0.98 

0.98 

0.98 

0.98 

0.97 

0.94 

1.03 

1.03 

1.02 

1 .01 

0.99 

0.96 

100 

0.92 

0.96 

0.98 

0.99 

0.99 

1.02 

1.17 

1.08 

1.07 

1 .04 

1 .02 

1.03 

300 

0.97 

0.97 

0.97 

0.97 

0.95 

- «- 

1.05 

1.04 

1 .04 

1.03 

1.02 

1 ,000 

0.94 

0.95 

0.95 

0.89 

0.94 

1.11 

1.07 

1.08 

1.12 

1 .10 

10,000 

0.50 

0.50 

0.49 

0.51 

0.79 

0.99 

1.12 

1 .10 

1.10 

1.30 

30,000 

1 .06 

1.57 

0.76 

0.69 

0.66 

0.71 

1.24 

0.95 

the  difference  Is  generally  within  10  percent  of  the  fission 
product  gamma  dose  except  for  the  yields  of  10MT  and  30MT.  For 
these  larger  weapons,  however,  the  overpressures  at  the  dose 
ranges  of  interest  are  generally  well  over  30  psi.  As  can  be  seen 
in  Table  1,  and  as  is  even  more  evident  from  listings  as  a function 
of  slant  range,  the  errors  arc  quite  systematic.  Thus,  If  desired, 
further  corrections  could  be  readily  developed  to  make  the  estimated 
error  still  closer  to  the  actual  error.  Such  corrections  would 
require  possibly  a 20  percent  to  50  percent  increase  in  calculation 
time  for  each  subroutine  call.  Use  of  this  multiple  approximation 
technique  is  not  untypical  of  this  approach,  where  the  error  bounds 
achieved  are  often  dependent  primarily  on  the  effort  expended  in 
developing  the  approximations. 
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APPENDIX  C 


FORTRAN  LISTING  OF  SUBROUTINE  WRADVN 

This  appendix  consists  of  the  entire  listing  of 
subroutine  WRADVN  of  the  DAMEVL  routine  of  the  TACWAR  code. 

Comment  cards  have  been  omitted,  and  more  detailed  explana- 
tions added  to  the  right  of  the  Fortran  statements  to  which 
they  apply. 

This  appendix  serves  three  purposes.  First,  it 
provides  an  example  of  the  degree  of  detail  required  for  a 
thorough  review  of  nuclear  effects  calculations.  Second, 
it  shows  the  implementation  of  reference  C--1  on  a modern 
digital  computer  and  explains  the  program  in  terms  under- 
standable to  readers  who  are  not  experts  in  Fortran  program- 
ming (although  some  knowledge  of  Fortran  is  necessary  to 
follow  the  program).  Third,  it  shows  the  implementation  of 
AP550  targeting  methodology  (calculation  of  weapon  radius 
only)  in  a more  advanced  programming  language  (Fortran  V) 
than  that  given  in  reference  C-2. 

i 
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Defense  Intelligence  Agency,  Physical  Vulnerability 
Handbook-Nuclear  Weapons  (U) , AP-550-1-2-69-INT , 

1 June  1969,  with  changes,  Confidential. 

Defense  Intelligence  Agency  Directorate  for  Intelligence, 
Mathematical  Background  and  Programming  Aids  for  the 
Physical  Vulnerability  System  for  Nuclear  Weapons, 
DI-550-27-74  , 1 November  1974  , with  Change  1. 
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APPENDIX  D 


TABULATED  COMPARISONS  OF 
SUBROUTINE  OFFCOV  WITH  EXACT 
NUMERICAL  INTEGRATION 


f 

J 

i 
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This  appendix  contains  the  results  of  comparing  the 
output  of  T AC WAR  subroutine  OFFCOV  with  numerical  integrations 
of  the  coverage  function.  Arguments  are  identified  as  follows: 
TAR  (target  radius/weapon  radius) , CEP  (circular  error  probable/ 
weapon  radius),  OFFSET  (aimpoint  offset/weapon  radius).  The 
ratio  of  the  coverages  (OFFCOV  algorithm/numerical  integration) 
is  also  tabulated. 
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Table  D- 3.  OFFCOV 
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Table  D-6 . OFFCOV 


1.75  0.096674  0.095103  0.983752  0.036944  0.035146  0.951328 

2.00  0.083465  0.082179  0-984590  0.035510  0.033784  0.951387 

2.50  0-058654  0.057879  0.986771  0.032291  0.030726  0.9S1541 

3.00  0.038096  0.037708  0.989820  0.028751  0.027362  0.951703 

4.00  0.012678  0.012670  0.999392  0.021393  0.020369  0.952155 
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